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ABSTRACT 

We report the most sensitive search yet made for solar-like oscillations. We observed the star 
a Cen A in B aimer- line equivalent widths over six nights with both the 3.9 m Anglo- Australian 
Telescope and the ESO 3.6 m telescope. We set an upper limit on oscillations of 1.4 times solar 
and found tentative evidence for p-mode oscillations. We also found a power excess at low 
frequencies which has the same slope and strength as the power seen from granulation in the 
Sun. We therefore suggest that we have made the first detection of temporal fluctuations due to 
granulation in a solar-like star. 

Key words: stars: oscillations - Sun: oscillations - Sun: atmosphere - stars: individual: 
a Cen A (HR 5459) - stars: individual: Procyon (HR 2943) 



1 INTRODUCTION 

Many attempts have been made to detect stellar analogues 
of the solar five-minute oscillations. As with helioseismol- 
ogy, it is hoped that the measurement of oscillation frequen- 
cies in other stars will place important constraints on stel- 
lar model parameters and provide a strong test of evolution- 
ary theory. However, despite several claims in the literature, 
it is fair to say that there has been no unambiguous detec- 
tion of solar-like oscillations in any star except the Sun (see 
reviews by Brown & Gilliland 1994; Kjeldsen & Bedding 
1995; Gautschy & Saio 1996; Heasley et al. 1996; Bedding 
& Kjeldsen 1998). 

Most searches for stellar oscillations have used one of 
two methods: (i) high-resolution spectroscopy, which aims 
to detect periodic Doppler shifts in spectral features, or 
(ii) differential CCD photometry, which has been used to 
search for fluctuations in the integrated luminosity of stars 
in the open cluster M67 (Gilliland et al. 1993). 

We have been using a new method to search for 
solar-like oscillations that involves measuring temperature 
changes via their effect on the equivalent widths (EWs) of 
the Balmer hydrogen lines. We found strong evidence for 
solar-like oscillations in the G subgiant r\ Boo (Kjeldsen et 
al. 1995; Bedding & Kjeldsen 1995), with frequency split- 
tings that were later found to agree with theoretical models 
(Christensen-Dalsgaard et al. 1995a,b; Guenther & Demar- 
que 1996). Since then, the improved luminosity estimate for 
?y Boo from Hipparcos measurements has given even better 
agreement (Bedding et al. 1998). However, a search for ve- 



locity oscillations in r\ Boo by Brown et al. (1997) failed to 
detect a signal, setting limits at a level below the value ex- 
pected on the basis of the Kjeldsen et al. result. More re- 
cently, Brown et al. (in preparation) have obtained a larger 
set of Doppler observations with lower noise which also fail 
to show convincing evidence for oscillations. 

Meanwhile, the equivalent-width method has now been 
used to detect oscillations in the spatially resolved Sun 
(Keller et al. 1998), in the 5 Scuti variable FG Vir (Viskum et 
al. 1998) and in the rapidly oscillating Ap star a Cir (Baldry 
et al. 1999). For the Sun it has been shown directly that the 
EW variations of Balmer lines arise because the line inten- 
sity oscillates much less than the continuum intensity (Ro- 
nanetal. 1991; Keller et al. 1998). 

We chose r\ Boo as the first solar-like target for the EW 
method because it was expected to have an oscillation am- 
plitude about five times greater than the Sun. This turned 
out to be the case (assuming the detection is real). We then 
turned to a Cen A, a more challenging target because of its 
smaller expected oscillation amplitude (comparable to solar; 
Bedding et al. 1996). Being a near twin of the Sun and ex- 
tremely nearby, this star is an obvious target to search for 
oscillations (e.g., Brown et al. 1994). Previous attempts us- 
ing Doppler methods were reviewed by Kjeldsen & Bedding 
(1995). They include two claimed detections at amplitudes 
4-6 times greater than solar (Gelly et al. 1986; Pottasch et 
al. 1992) and two negative results at amplitudes about 2- 
3 times solar (Brown & Gilliland 1990; Edmonds & Cram 
1995). 

Here, we report observations of a Cen A in Balmer-line 
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equivalent widths which set an upper limit of only 1 .4 times 
solar. We find tentative evidence forp-mode oscillations. We 
also find strong evidence for stellar granulation in the power 
spectrum, at a level consistent with that expected on the basis 
of solar observations. 



2 OBSERVATIONS 

We observed a Cen A over six nights in April 1995 from 
two sites: 

• At Siding Spring Observatory in Australia, we used the 
3.9-metre Anglo- Australian Telescope (AAT) with the coude 
echelle spectrograph (UCLES). We recorded three orders 
centred at Ha and three orders at H/3, which was possible 
thanks to the flexibility of the CCD controller. The weather 
was about 85% clear but transparency was variable. The ex- 
posure time was typically about 30 s, with a deadtime of 
23 s between exposures. This relatively long exposure time 
for such a bright target was achieved by spreading the light 
along the spectrograph slit by (i) defocussing the telescope 
and using a wide slit (5-10arcsec) and (ii) trailing the star 
backwards and forwards along the slit with a peak-to-peak 
amplitude of about 2arcsec and a period of about 4 s. The 
slit length was 14 arcsec. Each spectrum typically produced 
about 9.0 x 10 7 photons/A in the continuum near Hen and 
about 3.5 x 10 7 photons/A near H/3. 

• At the European Southern Observatory on La Silla 
in Chile, we used the ESO 3.6-metre telescope with the 
Cassegrain echelle spectrograph (CASPEC). We recorded 
three orders centred at Ha. The weather was 100% clear. 
The exposure time was typically about 10 s, with a deadtime 
of 12 s between exposures. The slit was 7 arcsec wide and 
10 arcsec long, with slightly less defocus than at the AAT 
and no trailing. Each spectrum typically produced about 
5.5 x 10 7 photons/A in the continuum near Ha. 

The B component of the a Cen system, which is 1.3 
magnitudes fainter than the A component, had a separation 
of 17.3 arcsec at the time of our observations, so there should 
not be any contamination. We nevertheless kept the spectro- 
graph slit aligned with the position angle of the binary sys- 
tem so that any light from the B component would be spa- 
tially separated on the detector. 

In addition to our primary target of a Cen A, we also 
observed Procyon for one hour at the start of each night, plus 
5 hours at both sites at the start of the fifth night. Finally, we 
observed the solar spectrum via the daytime sky (about 18 hr 
at AAT and 14 hr at ESO, spread over the six days). For both 
Procyon and the Sun, this amount of data turned out not to 
be sufficient to detect oscillations, but we did find evidence 
for granulation power (see Section 5). 

3 DATA PROCESSING 

Here we describe the steps involved in processing the data. 
Section 3.2 is specific to the equivalent-width method, while 



the other steps could apply, at least in part, to other types of 
observations (Doppler shift or photometry). 



3.1 Preliminary reduction 

(i) Correction of each CCD frame for bias by subtracting 
an average bias frame and then subtracting a constant that 
was measured from the overscan region of the CCD frame. 

(ii) Correction for CCD non-linearity. Measuring oscilla- 
tions at the ppm level requires that the detector be linear to 
the level of one part in 1000 or better. This cannot be taken 
for granted and our tests of different CCDs and controllers 
often reveal deviations from linearity of up to a few per cent. 
Unless correction is made for these effects, the extra noise 
will destroy any possibility of detecting oscillations. 

Both CCDs were Tektronix 1024 x 1024 chips, used at a con- 
version factor of about 12 electrons per ADU. The onset of 
saturation was about 450 000 electrons per pixel for the AAT 
and 250 000 for ESO. From our measurements, the linearity 
corrections for both were well approximated by the relation: 

ADU b s /ADUt r ue = 1 + a ADU bs 

-exp(ADU obs //3 + 7 ), (1) 

where ADU bs and ADUtruc are observed and true counts 
measured in ADU and (a, (3, 7) had values of (— 1 .63 x 10~ 7 , 
700, -59.33) for the AAT and (-2.35 x 10~ 7 , 700, -33.53) 
for ESO. The measurement error on this correction was one 
part in 3500. 

(iii) Correction for pixel-to-pixel variations in CCD sen- 
sitivity by dividing by an average flat-field exposure (all flat- 
field exposures were corrected for non-linearities before av- 
eraging). 

(iv) Subtraction of sky background, which was quite sub- 
stantial during twilight. The background was estimated from 
the regions at each end of the spectrograph slit, either side of 
the stellar spectrum in each echelle order. 

(v) Extraction of one-dimensional spectra. During this 
step, the seeing in each frame (i.e., the FWHM along the 
spectrograph slit) and the position of the star (i.e., light cen- 
troid along the slit) were recorded for use in decorrelation 
(see Sec. 3.4 below). 



3.2 Measuring equivalent widths 

Achieving high precision requires more than simply fitting a 
profile. The method described here was developed after try- 
ing several different approaches. By analogy with Str0mgren 
Hf3 photometry, we calculated the flux in three artificial fil- 
ters, one centred on the line (L) and the others on the con- 
tinuum both redward (R) and blueward (B) of the line. For 
each spectrum, the following steps were followed: 

(i) Placement of the three filters at their nominal wave- 
lengths. 
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(ii) Calculation of the three fluxes B, L and R* 

(iii) Adjustment of the slope of the spectrum so that R 
and B were equal. This was done by multiplying the spec- 
trum by a linear ramp. 

(iv) Re-calculation of the filter fluxes and hence of the 
equivalent width: W = (R - L) / R. 

(v) Repetition of steps (ii)-(iv) with the three filters at dif- 
ferent positions. Iteration to find the filter position that max- 
imized the value of W. The outputs were: W, the position 
of the line, the height of the continuum (from R), and slope 
of continuum. 

(vi) Repetition of steps (i)-(v) for four different filter 
widths. 

The result was four times series (Wi, W 2 , W 3 , W4), one for 
each filter width. 



3.3 Initial time series processing 

The quality of the data, as measured by the local scatter, var- 
ied considerably from hour to hour and night to night. The 
following procedure was applied to each night of data sepa- 
rately. 

(i) Clipping of each time series to remove outlying points 
(4ct clipping, where a is the local rms scatter). 

(ii) Calculation of weights for each time series. This in- 
volved assigning a weight to each data point that was in- 
versely proportional to the local variance (cr 2 ). 

(iii) Calculation of a w , the weighted rms scatter of each 
time series, using the weights just calculated. This was used 
to select the best filter width, i.e., the one which minimized 
g w . By using a weighted rms scatter, we do not give too 
much importance to the bad segments of the data. In practice, 
rather than choosing one filter width, we used a weighted 
combination. That is, we chose the powers a, b, c, d to mini- 
mize the weighted scatter on the time series W^W\W%Wf, 
where a + b + c + d= 1. 

This step and the subsequent ones rely on the fact that any 
oscillation signal will be much smaller than the rms scatter 
in the time series. Most of the scatter is due to noise and 
any method of reducing the scatter should be a good thing, 
although care must be taken not to destroy the signal or to 
introduce a spurious signal. 

3.4 Decorrelation of time series 

As well as measuring the parameter which is expected to 
contain the oscillation signal (W), we also monitored ex- 
tra parameters. The aim was to correct for instrumental and 
other non-stellar effects. For example, if W is correlated 
with the seeing, we would suspect some flaw in the reduc- 
tion procedure, since the stellar oscillation does not know 

* The flux in a filter is simply the total counts in the stellar spectrum after 
it has been multiplied by the filter function. 



what is happening in the Earth's atmosphere. By correlat- 
ing measured equivalent widths with seeing variations, one 
has a chance to remove the influence of seeing simply by 
subtracting that part of the signal which correlates with see- 
ing. This process of decorrelation, which can be repeated for 
other parameters (total light level, position on detector, slope 
of continuum, etc.), is very powerful but can also be quite 
dangerous if not done with care (see Gilliland et al. 1991 for 
a fuller discussion). 

Again, the process was applied to data from one night at 
a time. Performing decorrelation over shorter intervals runs 
the risk of moving power around and creating or destroying 
signal - simulations were useful to check these effects. 

Prior to decorrelation, we high-pass filtered the time se- 
ries to remove low-frequency variations. This was done by 
subtracting from the time series a smoothed version that was 
produced by convolution with a supergaussian envelope hav- 
ing a FWHM of 40 minutes. 

The AAT spectra included a strong Fe I line near H/3 at 
4383 A and we used the EW of this line as a decorrelation 
parameter. This line is temperature sensitive, with the oppo- 
site sign and similar amplitude to the Balmer lines (Bedding 
et al. 1996), so it should contain oscillation signal in anti- 
phase with the Balmer lines - this was allowed for in the 
analysis. We found that decorrelating against the Fe I EW 
gave a significant reduction in the scatter of both the Ha and 
H/3 EW time series. Unfortunately, the ESO spectra (around 
Ha) did not contain any strong iron lines. 

Figure 1 shows the final time series of EW measure- 
ments and weights. The data quality clearly varies signifi- 
cantly. Almost one third of the AAT spectra of a Cen A were 
discarded, although this is not quite as bad as it sounds be- 
cause many of these were taken well into morning twilight. 
Still, many night-time AAT spectra were given zero weight. 
The reason for these poor measurements is not clear, but it is 
probably related to seeing-induced variations in the point- 
spread-function through the wide slit (see below). Those 
AAT measurements that are useful have much lower scat- 
ter, and therefore higher weights, than ESO measurements 
because of the availability of the Fe I line for decorrelation. 



4 RESULTS AND SIMULATIONS 

The amplitude spectrum of the EW time series was calcu- 
lated as a weighted least-squares fit of sinusoids (Frandsen et 
al. 1995; Arentoft et al. 1998), scaled so that a sinusoid with 
an amplitude of 1 ppm produces a peak of 1 ppm. The power 
spectrum, the square of this amplitude spectrum, is shown in 
Figure 2 for the region in which oscillations are expected. No 
obvious excess of power is seen. We can fit a two-component 
noise model to the whole power spectrum which consists of 
(i) white (i.e., flat) noise at a level of 4.3 ppm in the am- 
plitude spectrum, and (ii) a non-white component, which is 
discussed in Section 5. 

The total noise at 2.3 mHz is 4.7 ppm in the amplitude 
spectrum, while that expected from photon statistics alone 
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Figure 1. The final time series of EW observations of a Cen A. The upper panel shows the square root of statistical weights (the inverse of the local rms 
scatter) in arbitrary units for all AAT measurements (both Ha and H/3 are shown). There are 8430 points, with 5773 having non-zero weight. For points with 
zero weight, we have added a random vertical scatter to make it easier to see their distibution in time. The middle panel shows the same for ESO data (Ha 
only); there are 8847 points, with 7999 having non-zero weight. The bottom panel shows all EW measurements having non-zero weight (13772 points). The 
peak-to-peak height of the individual modes we are trying to detect corresponds approximately to the size of the dots in the figure. 
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Figure 2. Power spectrum of equivalent-width observations of a Cen A in the region where signal would be expected. The right-hand plot shows the spectral 
window, with the same horizontal scale, calculated by taking the power spectrum of a sinusoid with an amplitude of 1 ppm and the same sampling and weights 
as the actual data. 



is 2.4 ppm. We are therefore a factor of 1.8 away from the 
photon noise, with the extra noise presumably coming from 
instrumental and/or atmospheric effects. The use of a wide 
slit, while advantageous in terms of exposure times and duty 
cycle, led to seeing-induced changes in the point-spread- 
function which may be the cause. 

4.1 Simulations 

We expect oscillations in a Cen A to produce a regular series 
of peaks in the power spectrum, with amplitudes modulated 
by a broad envelope centred at about 2.3 mHz (Kjeldsen & 



Bedding 1995). The average amplitude of oscillation modes 
near the peak of this envelope, as measured in Balmer line 
EW, is expected to be about 8 ppm, which is 1 .4 times that 
for the Sun (Bedding et al. 1996). To investigate whether 
an oscillation signal may be present in our data, we have 
generated simulated time series consisting of artificial sig- 
nal plus noise. Each simulated series had exactly the same 
sampling function and allocation of statistical weights as the 
real data. The injected signal contained sinusoids at the fre- 
quencies calculated by Edmonds et al. ( 1 992) for modes with 
I = to 3, modulated by a broad solar-like envelope cen- 
tred at 2.3 mHz and with a central height of 8 ppm. In each 
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simulation, the phases of the oscillation modes were cho- 
sen at random and the amplitudes were randomized about 
their average values. All these characteristics were chosen to 
imitate as closely as possible the stochastic nature of oscilla- 
tions in the Sun. Before calculating each power spectrum, we 
added normally-distributed noise to the time series, so as to 
produce a noise level in the amplitude spectrum of 4.7 ppm 
(consistent with the actual data). 

Some results are shown in Figure 3. The top panel 
shows a simulation without any added noise; the randomiza- 
tion of mode amplitudes within the broad envelope is clear. 
The next five panels show simulations with noise included 
(each with different randomization of noise, mode ampli- 
tudes and mode phases), while the bottom panel shows the 
actual data. 

It is interesting to note that some of the signal peaks 
in the simulations have been strengthened significantly by 
constructive interference with noise peaks. For example, a 
signal peak of 8 ppm which happens to be in phase with a 
2<T noise peak (2 x 4.7 ppm) will produce a peak in power 
of 300 ppm 2 . This illustrates the point made by Kjeldsen & 
Bedding (1995; Appendix A.2): the effects of noise must be 
taken into account when estimating the amplitude of a sig- 
nal. 

In these five simulations we can see that a signal of 
8 ppm is sometimes obvious but sometimes not. To quan- 
tify this, we have looked at the eight strongest peaks in each 
power spectrum in the frequency range 1600-3000 ^Hz. 
These are marked by vertical lines in Figure 3. Dotted lines 
show peaks that coincide with input frequencies or their 
1/day aliases to within ±1.3 fiHz (which is twice the rms 
scatter on the differences between input and measured fre- 
quencies). The remaining peaks, marked by dashed lines, are 
assumed to be due to noise. The results show that 2-5 of the 
eight strongest peaks coincided with input frequencies and 
a further 0-2 were 1/day aliases (the mean values for these 
were 4.0 and 1.0, respectively). 

The same information is shown in Figure 4 in the form 
of echelle diagrams. The top panel shows all input frequen- 
cies from Edmonds et al. (1992), although most of these are 
given very small amplitudes in the simulations. The symbols 
indicate different I values: (squares), 1 (triangles), 2 (dia- 
monds) and 3 (crosses). The next five panels show the eight 
strongest peaks in each simulation. The solid lines show the 
loci of input model frequencies. Each peak in Figure 3 is 
identified either as noise (asterisks) or as one of the input fre- 
quencies. In the latter case, some peaks are shifted by 1/day 
(11.57^Hz), as indicated by dotted lines. 



3000 
2500 

2000 

1500 
3000 

2500 

2000 

1500 
3000 

2500 
2000 

1500 
U 3000 

a: 
a. 

2500 

o 
c 

3 2000 
cr 

£ 1500 
3000 

2500 

2000 

1500 
3000 

2500 

2000 

1500 
3000 

2500 

2000 

1500 





-20 











20 40 60 

Frequency modulo b±v (/iHz) 



B0 



Figure 4. Echelle diagrams based on Figure 3. The top panel shows the 
input model frequencies, plotted using Av = 108.0 /iHz. The symbols 
indicate different £ values: (squares), 1 (triangles), 2 (diamonds) and 3 
(crosses). The next five panels show the eight strongest peaks in each simu- 
lation. The solid lines show the locus of input model frequencies. Each peak 
is identified either as noise (asterisks) or as one of the input frequencies. In 
the latter case, some peaks are shifted by 1/day (11.57 /iHz), as indicated by 
dotted lines. The bottom panel shows the eight strongest peaks in the actual 
data, with Av = 107.0 /iHz. 



4.2 Mode identifications 

Given that some of the eight strongest peaks in our power 
spectrum of a Cen A may be real, we can attempt to identify 
the modes to which they correspond. Mode frequencies for 
low-degree and high-order oscillations in the Sun and other 
solar-like stars are well approximated by the asymptotic re- 
lation: 



v{n,t) = A*/(n+ \l + e) - £(£ + 1)D . (2) 

Here n and I are integers which define the radial order and 
angular degree of the mode, respectively; /Av (the so-called 
large separation) reflects the average stellar density, D is 
sensitive to the sound speed near the core and e is sensitive 
to the surface layers. Values for these three parameters in the 
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Figure 3. Simulated power spectra for a Cen A, using the same sampling times and data weights as the actual observations. The top panel shows a simulation 
without noise, while the next five panels show simulations with noise. The eight strongest peaks are marked by vertical lines, with dotted lines indicating those 
coinciding with input frequencies or their 1/day aliases, and dashed lines showing noise peaks (see also Figure 4). The bottom panel shows the actual data 
from a Cen A, with the eight strongest peaks being marked (see Table 1). 



Sun are: 

Av = 135.12 ± 0.18 [iHz 

D = 1.50 ± 0.03 fjHz 

e = 1.46 ±0.03 

We obtained these values by fitting Equation 2 to solar fre- 
quency measurements for n = 17-25 and I = 0-2 by 
Frohlich et al. (1997). 



The curvature in the Edmonds et al. (1992) frequencies 
for a Cen A, obvious in Figure 4, indicates a departure from 
the asymptotic theory. Curvature is also predicted by models 
of the Sun, at a level larger than actually observed, reflect- 
ing the difficulty in modelling the solar surface. Therefore, 
the real oscillation frequencies in a Cen A are also likely 
to show less curvature than the model by Edmonds et al. 
This does not affect the conclusions of the simulations, but 
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Table 1. Frequencies of the eight strongest peaks in the a Cen power spectrum and possible mode identifications (n, £). 



Frequency Possible identifications 



(/iHz) 


Case 1 


Case 2 


Case 3 


Case 4 


1687.79 


i/(13,3) - 1/d? 


noise 


"/(15,1) 


1/(14,2) 


2145.76 


i/(18, 1) 


i/(18,2) 


j/(20,0) 


1/(19,1) 


2192.55 


i/(18,2) 


i/(18,3) 


1/(20,1) 


1/(19,2) 


2467.68 


1/(21, 1) 


i/(21,2) 


noise 


noise 


25f2.35 


i/(21,2) 


i/(21,3) 


noise 


i/(23, 0) + 1/d? 


2648.48 


i/(22,3) - 1/d? 


noise 


i/(25,0) 


1/(24, 1) 


2799.92 


i/(24, 1) + 1/d? 


i/(24,2) + 1/d? 


noise 


noise 


2896.28 


i/(25, 1) 


f(25,2) 


1/(27, 1) 


1/(26,2) 




106.94 ± 0.13 


106.99 ±0.16 


100.77 ±0.17 


100.77 ±0.17 


D //uHz 


2.05 ± 0.05 


1.35 ±0.09 


1.95 ±0.29 


1.07 ±0.09 


e 


1.61 ± 0.03 


1.14 ±0.04 


1.29 ±0.04 


1.81 ±0.04 



is important when we try to assign modes to our observed 
frequencies. 

The bottom panel of Figure 3 shows the eight strongest 
peaks in the actual data. Assuming our representation of the 
mode amplitudes in a Cen A is realistic, we would expect 
some of these to correspond to actual frequencies. We have 
therefore searched for values of Av which fit the observed 
frequencies, as follows. We calculated all pairwise differ- 
ences and looked for a value of Av which agreed with as 
many of these as possible (within ±0.8 /iHz). The only two 
values of Av in the range 96-1 16 /iHz which gave four or 
more coincidences were 100.8 /iHz (6 pairs) and 107.0 /iHz 
(5 pairs). We also checked how well randomly chosen fre- 
quencies gave similar coincidences. On the basis of 25 sim- 
ulations, we found that 6 coincidences occurred only once. 

For both of the values of Av found above, we attempted 
to identify each of the eight observed frequencies by fitting 
to Equation 2. There are two possible identifications for each 
value of Av that have Da in a physically realistic range (0.7- 
3.0 /iHz; Christensen-Dalsgaard 1993) and these are given in 
Table 1. We chose n such that e lay in the range 1-2. 

We stress that we do not claim to have detected oscil- 
lations but that, if we have, one of the four cases in the ta- 
ble represents the most likely description. In that case, the 
implied amplitude is 7-8 ppm (1.3 times solar). If none of 
the eight frequencies is real, the oscillations must be below 
7 ppm (1.2 times solar). Either way, these observations repre- 
sent the most sensitive search yet made for stellar variability. 



5 GRANULATION POWER 

The power spectrum of full-disk intensity measurements of 
the Sun shows a sloping background due to granulation and 
other surface structure (e.g., Rabello Soares et al. 1997). This 
power background reflects surface temperature fluctuations, 
so we expect it also to be present in EW measurements. In- 
deed, our power spectrum of a Cen A (Figure 5) shows an 
excess at low frequencies. The same is true for our observa- 



tions of Procyon and the Sun, although there are much less 
data. 

Note that the power density scales on the right-hand 
axes were obtained by multiplying the left-hand scales (de- 
fined in Section 4) by the effective timespan of the observa- 
tions (see Appendix A.l of Kjeldsen & Bedding 1995). For 
a Cen A, for example, the observations spanned 145 hr from 
start to finish. However, the measurements were given differ- 
ent weights, so we calculated the effective timespan by inte- 
grating under the weighted spectral window (right panel in 
Figure 2); the result was 85.5 hr. It is important to note that 
published power density spectra are often calculated using 
different versions of Parseval's theorem and must be multi- 
plied by either two or four to be compared with our definition 
of power density. 

To investigate whether the low-frequency power excess 
observed in a Cen A could arise from stellar granulation, 
we must make two corrections. Firstly, we must subtract the 
contribution from white noise (photon noise, instrumental 
and atmospheric effects), which we measure to have a level 
of 7.3 ± 0.2 ppm 2 //iHz. This is shown by the horizontal line 
in Figure 5. Secondly, at frequencies below about 800 /iHz 
there is a deficit of power due to the high-pass filter (see 
Section 3.4). 

The thick solid line in Figure 6 shows the power density 
from our EW measurements of a Cen A after subtraction of 
the constant white noise term, correction for high-pass filter- 
ing and smoothing. The two thinner solid lines on either side 
show the ±2cr errors in the two corrections. 

The power is quite uncertain at low frequencies due to 
uncertainties in the removal of the effects of the high-pass 
filter, and at high frequencies due to uncertainties in the sub- 
traction of white noise component. At intermediate frequen- 
cies (600-3000 /iHz), our measurement should be accurate 
and in this regime the power is linear and well described by 
the following relation: 



log (power density) = G + P log 
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Figure 5. Power spectra of equivalent-width observations of a Cen A, Procyon and the Sun. In each case, the mean level of white noise is shown by the 
horizontal white line. 
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Figure 6. Comparison of power density spectra. The thick solid curve shows the observed EW power density from our AAT7ESO observations of a Cen A, 
after subtraction of white noise and after correction at low frequencies for the effect of the high-pass filter. The two thinner solid lines on either side show 
the ±2cr errors in these two corrections. Dashed curve: same as thick solid curve, but without correction for the high-pass filter. The symbols show EW 
measurements from our AAT/ESO observations for Procyon and the Sun, with 2<r error bars. The dotted line shows our estimate of the EW granulation power 
density in the Sun, scaled from published intensity measurements by the VIRGO LOI instrument (see text). 
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Figure 7. Example time series for EW measurements of a Cen A based of the non-white component of the power spectrum. 



where 

G = 0.85 ±0.09 
P = -1.46 ±0.15 

and the power density is measured in ppm 2 //iHz. The slope 
of the background power is very similar to that produced by 
solar granulation (e.g., Rabello Soares et al. 1997). To com- 
pare the amount of power in our EW measurements with 
published measurements of the total solar intensity, we must 
estimate the conversion factor between the two observing 
methods. For the Sun, a fractional change in intensity of 
1 ppm that is caused by a temperature change corresponds 
to a change in EW of the Balmer lines of 1 .5 ppm (Bedding 
et al. 1996). We assume the same conversion factor applies 
to temperature fluctuations from granulation. Furthermore, 
the amount of limb darkening in the EW of Balmer lines 
- and hence the signal from granulation - is greater than 
in total intensity (Bedding et al. 1996) by a factor which 
we estimate to be 1.15. Hence, published measurements of 
the power density of solar intensity should be multiplied by 
(1.5) 2 (1.15) 2 = 3.0. 

We have applied this conversion to the recent mea- 
surement of full-disk solar intensity measurements made 
with the LOI (Luminosity Oscillations Imager) part of the 
VIRGO instrument on the SOHO spacecraft (Fig. 2 of Ap- 
pourchaux et al. 1997). The result is shown as the dotted line 
in Figure 6. Note that it was also necessary to multiply the 
LOI values by four to bring them into line with our defini- 
tion of power density (T. Appourchaux, private communi- 
cation). The agreement with our observations of a Cen A is 
excellent, giving strong evidence that we have detected stel- 
lar granulation. 

Our observations of Procyon and the daytime sky also 
produced power spectra with excesses at low frequencies 
(Figure 5). These power levels are shown in Figure 6, with 
2cr error bars. Our solar measurement is in reasonable agree- 
ment with the LOI data, although slightly too low. The gran- 
ulation power in Procyon appears to be greater than that in 
a Cen A by a factor of 2.0, although this result is only at the 
2cr level. 

To illustrate the size of the granulation signal from 
a Cen A, we show in Figure 7 an example time series, sam- 



pled at 50 s intervals, based on the non-white component of 
the power spectrum (the thick line in Figure 6). This is how 
the EW signal from a Cen would look without photon noise 
or measurement errors. 

When one is trying to detect oscillations, these fluctu- 
ations in EW due to stellar granulation represent a funda- 
mental noise source. For example, suppose that the weather 
had been 100% clear and that all of our EW measurements 
(Ha and H/3) were photon-noise limited. Then we would 
have expected a white noise level in the amplitude spectrum 
(from photon noise) of 1 .4 ppm and a noise level from gran- 
ulation at 2.3 mHz of 1.9 ppm. The total noise at 2.3 mHz 
would then have been 2.4 ppm, about half the value that we 
actually measured. In other words, we were a factor of two 
away from the limit set by photon and granulation noise, pre- 
sumably due to instrumental and atmospheric effects. Had 
we achieved this limit, a signal with solar amplitude (about 
6 ppm) would have been detected. 

Given a CCD with sufficiently fast readout, it would be 
possible to record a much larger part of the spectrum and 
hence to measure EWs of many more lines at a higher duty 
cycle (Bedding et al. 1996). The photon noise would then be 
reduced by up to a factor of about two, to which granulation 
noise must still be added. This represents the fundamental 
limit for the EW method. 

Granulation is also a fundamental limit for intensity and 
colour measurements, but much less critically for Doppler 
shift measurements, since the background power in the solar 
velocity spectrum is a tiny fraction of the oscillation ampli- 
tudes (Gabriel et al. 1997). 



6 CONCLUSIONS 

Our observations of a Cen A represent the most sensitive 
search yet made for solar-like oscillations. We can set a firm 
upper limit on the amplitudes of 1 .4 times solar, which is ap- 
proximately the level at which oscillations are expected. We 
find tentative evidence for p modes in the form of a regular 
series of peaks. If these do not correspond to real oscillations 
then the upper limit on amplitude becomes 1 .2 times solar. 
We find an excess of power at low frequencies in 
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a Cen A which has the same slope and strength as power 
from granulation in the Sun. We therefore suggest that we 
have made the first detection of granulation power in a solar- 
like star. 
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